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We perform a quenched computation of the glue momentum fraction in the pion. Different dis- 
cretizations of the gluonic energy-momentum tensor are studied on the lattice for that purpose. 
We discuss some implications based on the momentum sum rule. Finally we point out promising 
applications of the techniques developed here. 
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1. Introduction 

A striking feature of QCD is the large contribution of gluons to the mass and momentum of 
hadrons, so it is of fundamental interest to calculate the contributions of gluons from first principles 
using lattice QCD. The first moments 

Wf(^') = Lf^uAs!o^dx{f{x,q')+f{x,q^)] (1.1) 
Wg(^^) = SQxdxg{x,q^) (1.2) 

of the quark and gluon distribution functions /(x), f{x) (f = u,d,s,.. .) and g{x) acquire a precise 
field-theoretic meaning via the operator product expansion in QCD. They satisfy the well-known 
momentum sum rule (MSR) {x)f{q^) + {x)g{q^) = 1 and are related to the corresponding contribu- 
tions to the energy-momentum tensor T^y evaluated on the hadronic state. Separating the traceless 
part T^iy from the trace part S for gluons, denoted 'g', and quarks, denoted 'f, T^y has the explicit 
form 

T^,y = Tl^ + T\^ + \5^y{S^ + S'), (1.3) 

— 1 f\ F'^ F'^ — F" F" n 4^ 

^ jiv pa^ pa ^iia^vai y^-^/ 

T^v = iLfWfD^ YvVf + ¥fDy Y^Yf - ^5^v¥fDp YpYf, (1-5) 

= p{g)/{2g) f^^f;;^, s' = [\ + YMmfVm (i-6) 

where D^=D^ -D^,^{g) = -b^g^ + ... is the beta-function with bo = {^N - jNf) (471)"^, Ym{g) 
is the anomalous dimension of the mass operator, and all expressions are written in Euclidean 
space. The gauge action in this notation is ^F^yF^y. For an on-shell particle with four-momentum 
p = {iEp,p), Ep = + p^, we have the relations 

{^,p\Jd'zTll{z) = [£p - |MV£p] (x)f,g, (1.7) 
{^,p\Jd'zS''^{z) \^,p) = (M^/Ep) bi^„ (1.8) 
{x)f+{x)g = bf + bg = \, (1.9) 

where states are normalized according to (pip) = 1. The renormalization of (jc)f,g will be discussed 
in section ^ 

Equation shows that in the infinite momentum frame, where Ep ^ P ^ oo, l^x)g represents 
the momentum fraction arising from gluons, and calculating {x)g is the main goal of this work. In 
the rest frame, the gluon contribution of Eq. 1.7 to the hadron mass is ^M{x)g From Eq. 1.8 



in the rest frame, the contribution of the trace anomaly to the hadron mass is \bgM [|l|], and we 
have performed the first step to calculate this matrix element as well Finally, we remark that in 
thermodynamics at temperature T the energy density e and pressure P are given by 

e-3P={S)T-{S)o, e + P=l{Too)T. (1-10) 

The pioneering calculation of the glue momentum fraction (in the proton) was presented at the 
LATTICE96 conference In the present study we treat the case of "heavy pions" with masses 
in the range 600 MeV < Mj^ < 1060MeV. Our improved techniques, applied here in the quenched 
approximation, are applicable in full QCD calculations, and to the case of the proton. 
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2. Discretization and normalization of the gluonic operators 

We use the Wilson gluon action -TLr,/j/vTr{l —P^y{x)} and the Wilson fermion action 
There are two distinct ways ||^ to discretize the Euclidean gluonic operator Tqq = ^(— E*^' • W + 
B" • B^') and the trace anomaly = ^ (E" • E" + B" • B") on a hypercubic lattice. 

The first, denoted 'bp' for bare-plaquette, uses a sum of bare plaquettes P^y around a body- 
centered point xq = X + \a'^^;^ fx: 



ro[;(x0) EE ;^''P(go)-^£^oo(0^,l^) L Tr{l-Po,^i,^(x + aA2;;+aw3^)} (2.1) 

°" ^0 ^ co,A=0,l 



5''P(xo) ^XvP(5o);^^I I Tr{l-Po,^i,^(x + aA2^ + a«3^)} (2.2) 
fliogasa ^(B,A=o,i 

where /i^ is the image under permutation ^ of jJ. and soq{ij.,v) is 1 if (/i, v) are both spatial and 
-1 otherwise. Other diagonal elements of the energy-momentum tensor are obtained by letting the 
other directions play the role of time. When summed over a time-slice, these expressions simplify 
to Eq. 7 of Ref. [||]. The other form ('bare clover') is 



Toq{x) = 2 ^^Tr l^iFok) -l^{Fti) 

So ^ k k<i ■ 

dloga 



(2.3) 
(2.4) 



where Fnv{x) is the clover-shaped discretization of the field-strength tensor (see [§]). This form 
allows for the discretizations of off-diagonal elements of T^y as well. Each of the normalization 



factors Zg(go), X (go) and zrfeo) in Eq. (|11|2J) is of the form 1 +0(g^). The factor Zg(go) is 



dictated by an exact lattice sum-rule for the Wilson gauge action and is known with a precision of 



about 1% (see [ ]10| ] and Refs. therein). To obtain the absolute normalization of other discretizations, 
it is sufficient to compute their normalization x{go) relative to that of the bare plaquette. 

New versions of the gluonic operators are obtained by replacing the Unk variables in Eq. (p|- 



2.4| ) by 'smeared' versions thereof. The HYP form of smearing ||9|] is particularly attractive in the 
present context in that it preserves the symmetry among all four Euclidean directions. 

Our criteria for the choice of the discretization are to maximize the signal-to-noise ratio, mini- 
mize cutoff effects, and preserve locality as much as possible. We studied the signal-to-noise 
ratio for four different discretizations by comparing the variance of the entropy at temperature 
T = 1/Lo = l.llTc 0, which is proportional to the expectation value of T^x^ooix), on Lq x 
lattices of fixed physical size. 

The resulting variances for the plaquette and clover discretizations with bare and HYP links 
are shown in Table [T[ We find dramatic differences between the discretizations, with the HYP- 
clover operator reducing the variance by almost two orders of magnitude as compared to the bare 
plaquette operator. Variance reduction comes at the cost of a certain loss of locality, since the HYP 
plaquette and HYP-clover operators have extent 3a and 4a respectively. As Fig. (|l|, left panel) 
illustrates, the variance of ^oo(-^) grows like a^'*, but the prefactor is non-universal and depends 
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relative variance 


normalization 


bare HYP 


bare HYP 


j3 = 6.000 

6x 16^ 


plaq. 

clover 


26.42(71) 0.6518(43) 
3.85(11) 0.3049(41) 


1.000 0.5489(68) 
2.184(67) 0.613(20) 


j3 = 6.093 

7 x203 


plaq. 
clover 


45.5(2.5) 0.991(11) 
6.20(14) 0.4393(40) 


1.000 0.546(14) 
2.026(52) 0.585(16) 


/3 =6.180 
8 x223 


plaq. 
clover 


81.5(5.8) 1.940(26) 
11.84(33) 0.8280(85) 


1.000 0.596(20) 
2.092(68) 0.636(22) 


j3 = 6.295 

9x24^ 


plaq. 
clover 


157(16) 3.080(50) 
19.74(75) 1.268(16) 


1.000 0.563(28) 
1.928(93) 0.598(30) 


j3 = 6.408 
11 x293 


plaq. 
clover 


346(81) 7.67(29) 
43.3(3.6) 3.042(78) 


1.000 0.617(68) 
1.91(20) 0.654(75) 



Table 1: Left : the relative variance, {ff^) / {&)^ — 1, of the operators & = Y.^ T^o{^) on lattices of matched 
physical sizes for different discretizations described in the text. Right: the normalization (go, fl/i^o) of the 
operator relative to the bare plaquette, determined on the same lattice. 



Wilson action P=6.0 K=0.1515 




Figure 1: Left : the relative variance, - 1, of the operator = 'LxTIq{x) as Hsted in Table 1. 

The lines of the form y — ax^ are drawn to guide the eye. Right: The effective gluonic momentum fraction 

(it) 

in a heavy pion, M^r ~ 1060MeV, using the HYP-plaquette discretization. 



strongly on the discretization. As a compromise between locality and variance reduction, from now 
on we work with the HYP-plaquette operator. In we performed a check of its discretization 
errors by computing the dependence of ;^ on a/Lo at j3 = 6. We showed that the dependence of x 
ona/Lo is mild and statistically consistent with zero for Lo/a>6. 

The X normalization factors are given in Tab. |l| for the four discretizations. We have thus deter- 
mined the non-perturbative normalization of Tqq in a range of lattice spacings [||], 5 < ro/a < 10. 
As a check of the correct normalization of the chosen HYP-plaquette operator, we computed its 

expectation value on the lightest scalar glueball. In that case, we know that the momentum fraction 

(c) 

carried by the glue is one, and indeed we find {xYg'' = 1.16(18). 
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Mn (Mev) 


32-123 


32-163 


48-163 


24^ 


1060(10) 


0.39(6)23091 


0.29(9)7113 


0.40(8)8331 


0.34(9)1048 


891(9) 








0.36(8)3066 


624(6) 








0.58(16)2538 



( 71) 

Table 2: The glue momentum fraction (jc)g ' in the pion. The integer in each subscript denotes the number 
of configurations used. 



3. The glue momentum fraction in the pion 



We consider a triplet of Wilson quarks, labeled u,d,s, with periodic boundary conditions in 
all directions and with common K" = 0.1515, 0.1530 and 0.1550 corresponding to pion masses 
approximately 1060, 890 and 620 MeV on the jS = 6 lattices 32 - 123, 32 . 153^ 4g . ^53 24'^. We 
use ro = 0.5fm to set the scale. To calculate the gluonic momentum fraction in the pion, we define 
the effective momentum fraction 



/g.eff 



3M;, lAol 



x;j:oeAo 



Iy(7(0)r;;f;(xo)j(^,y)) ,^hp 



;¥,y')> 



^00 fe)) 



(3.1) 



.mm 
' 



Lo 
' 2 



mm_ Lo _|_ min 



. ,Lo-xr-a}. For large Lo and.^^' ^ i^^ 



where Ao = {x\ 

As a source field for the pion, we use the isovector pseudoscalar density j{x) = d{x)Y5u{x). Its two- 
point function is positive on every configuration, for each of which we do 12 inversions correspon- 
ding to Dirac and color indices. On a 24^* lattice, we take advantage of the symmetry between all 
directions to perform these inversions at the points ^(6,6,6,6) for ^ = 0, 1,2,3 and symmetrize 



expression ( [3.1[ ) with respect to all directions, so that ^^^^ T^^{x) vanishes on every configuration. 



The right panel of Fig. |T] shows our stable plateaus for (jc)g jff at large values of .x™" for three lattice 
sizes, and the results are summarized in Tab. ^ 



The lattice equivalent of Eq. |L9| has been derived for Wilson lattice QCD in []11[]. It reads 



1 + 0{a') = (x)g + {x)f, {x)f = Zf(go)(x)^^ (x)g = Zg(go)(x 



\bare 



(3.2) 



The particular discretization of Tqq appearing in the lattice sum rule is the bare plaquette one. The 
normalization factor Zg(go) and its fermionic analog Zf(go) are given by the so-called anisotropy 
coefficients used in thermodynamic studies. The MSR thus holds even in an 'unrenormalized form' 
where (A')f g do not separately have a continuum limit; their renormalization is performed in the next 
section. 



4. Renormalization of {x)c 



Recall that, in QCD, the renormalization pattern in the singlet sector reads [|T2p 



>oo(Ai)" 
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1-Z 



gg 



1-Zff 
Zff 



T'ooteo) 



(4.1) 
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provided (go) are normalized so that Eqs. ([P7|-|r^) hold. In lattice regularization, this requires 
the scheme-independent Zg(go) and Zf (go) factors, whereas Zgg and Zff are scheme-dependent func- 
tions of (a/ijgo)- The renormalization group equation then takes the form 



(^)g(M') 
(^)f(M') 



-(M) 



Cgg(f) 
-Cgg{g) 



-Cff(g) 
Cff(g) 



(^)g(At^ 
(x)f(At' 



(4.2) 



JVf_ 
12;r2 • 



respectively [13]. Be- 



with nd^ log[Zgg + Zff - 1] = -g2 [cgg + Cff] and Cgg^ff (g = 0) 

sides the zero-mode Too, the linear combination [1 + x{pL)]T'^{}i.) + T(/x)roo(/x) renormalizes 
multiplicatively with anomalous dimension — g^[cff + Cgg], where jxd^T = — g^[(cff + Cgg)T + Cff]. 
Note that the asymptotic glue momentum fraction is given by Cff(0)/[cff(0) +Cgg(0)] = Zgg(oo) = 
l-Zff(oo) = -T(oo) = 16/[16 + 3M-]. 

In the quenched approximation, Zgg = 1 due to the absence of quark loops [[l^, |l6|]. This 
implies that the singlet quark operator renormalizes multiplicatively and with the same Z-factor, 



computed non-perturbatively in [14], as the non-singlet operator: 



(^)g(M^) = (^)g + [l-Zff(aAi,go)] (x)f, 
(x)f(At^) = Zs{aix,go) (x)f. 



(4.3) 
(4.4) 



Disregarding disconnected diagrams, {x)^^'^ has been computed in [15] at the bare parameters 
(j8 = 6, K- = 0.1530). The product Zff(aAt,go)Zf(go) = 0.99(4) (in the MS-scheme at /i = 2GeV) 
is known from [[T^, |T5|]. The factor Zf(go) = 1 + 0(go) is as yet unknown beyond tree level. If we 
allow for a conservative error, based on the typical size of one-loop corrections, Zf(go) = 1.0(2), 
our final result is 



(;c)J^'(At^ = 4GeV2) =0.37(8)(12) {M^ = 890MeV), 



(4.5) 



where the first error is statistical and the second comes from the uncertainty in Zf(go). We observe 
no significant quark-mass dependence, and given the uncertainties, we do not attempt a chiral 
extrapolation. 

Our result for the glue momentum fraction in a (heavy) pion is compatible with phenomeno- 
logical determinations [jTsj, [I^, {xYf = 0.38(5) at = 4GeV2, based on Drell-Yan, prompt 
photoproduction, and the model assumption that sea quarks carry 10-20% of the momentum. The 
agreement suggests a mild quark-mass dependence, but only a calculation in full QCD and at 
smaller masses can substantiate this. Our result at = 4GeV^ lies clearly below the Ni = 3 
asymptotic glue momentum fraction of 0.64. The fact that our result and the valence quark mo- 



mentum fraction, computed in [15], add up to 0.99(8)(12) suggests that the omitted disconnected 
diagrams are small. 



5. Conclusion 

We have computed the glue momentum fraction {x)g in a pion of mass 0.6GeV < Mji < 
1.06GeV using quenched lattice QCD simulations. We find 37(8)(I2)% at = 2GeV, a resuh 
compatible with phenomenological determinations [ |T8| , p^ ]. 
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Although it appears difficult to achieve precision at the percent level, the present method is 
applicable to full QCD with dynamical quarks. Presently the larger uncertainty comes from the 
normalization of the quark contribution to the renormalized (x)g, and could be reduced significantly 
by a one-loop calculation. 

We also evaluated the bare trace anomaly contribution to the pion's mass in the same frame- 
work The counterterm remains to be calculated, but it will ultimately be preferable to use chiral 
fermions to avoid mixing with the lower dimensional fermion operator. 

Finally, we remark that the freedom of choosing a numerically advantageous discretization 
of T has not been fully exploited in previous lattice simulations. The improvement that was 
essential in the present computation of the pion momentum fraction can be carried over to fully 
dynamical calculations and the exploration of other observables, such as the gluon contribution to 
the nucleon spin. It is also particularly promising for thermodynamic studies of pressure, energy 
density and transport coefficients. 
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